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Abstract 

We study a model of mass-bearing coagulating planar Brownian particles. Coagu- 
lation occurs when two particles are within a distance of order e. We assume that the 
initial number of particles N is of order | log er| . Under suitable assumptions of the 
initial distribution of particles and the microscopic coagulation propensities, we show 
that the macroscopic particle densities satisfy a Smoluchowski-type equation. 

1 Introduction 

A colloid consists of a large number of small particles that are suspended in an environment 
of far smaller and more numerous molecules. Large numbers of molecules bombard each 
particle, and random fluctuations among these collisions tend to give rise to an Ornstein- 
Uhlenbeck motion of the particle, in which its velocity is forced by a Brownian motion, 
with a drag force acting in the direction opposite to its velocity. On a long time scale, 
the colloidal particles move according to Brownian motions, because an Ornstein-Uhlenbeck 
process approximates such a motion over a long period of time. The particles of a colloid 
may also be liable to interact. In [2], we studied a model of a colloid in which this means of 
interaction took the form of a coagulation, this reaction being liable to take place between 
a pair of particles if they come to lie close enough to one another. The density of particles 
at the initial time was chosen so that the dynamics occur in a regime of mean free path, 
wherein a typical particle meets a bounded number of other particles in a unit of time. 
Speaking in rough terms, this choice of scaling causes the effects of diffusion and interaction 
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on the macroscopic evolution of the system to be comparable. In common with much of 
non-equilibrium statistical mechanics, we interpret the macroscopic behaviour of the system 
in terms of the evolution of a small number of thermodynamic parameters, in this case, the 
density of particles of a given mass, as a function of macroscopic space and time. In [2], we 
proved that, when the initial number of particles is chosen to be high, this density typically 
evolves as the solution of the Smoluchowski system of PDE, 

^(s, t) = d{n) Af n (x, t) + Qi(f)(x, t) - Q n 2 (f)(x, t). n = 1, 2, . . . (1.1) 

The first term on the right-hand-side of (1.1) corresponds to the diffusion among particles of 
mass n, with d(n) being one-half of the diffusion rate of such particles. The terms in (1.1) 
corresponding to the interaction of pairs of particles are given by the gain term 

n 

QUf) (x, t) = P( m > n - m)f m (x, t)f n - m (x, t), (1.2) 

m=l 

and the loss term 

oo 

Q n 2 {f) = 2f n (x, t) Pirn, n)f m (x, t). (1.3) 

m=l 

Here, the collection of constants (3 : N 2 — ► (0, oo) quantify the macroscopic propensity of 
mass at a pair of values to combine. 

The arguments of |5j were valid in the case where the dimension d of the system was 
assumed to be at least three. The question of the behaviour of such a system in two dimen- 
sions is significantly different, and it is this topic that we address in this paper. We now 
turn to describe the model in more detail, after which, we will discuss the ways in which the 
two-dimensional case differs from that of higher dimensions. 

We will be working with a collection of microscopic models, each model carrying an index 
iVeff, this being the total number of particles present in the system at the initial time. 
Each of these N particles is independently assigned a random integer mass and placed at the 
initial time at a random location whose law depends on that mass. More precisely, we will be 
describing the state of the system at any given moment in time by a configuration, by which 
we mean a map q : I q — > IR 2 x N, whose domain I q is some finite set of a countable index 
set I. That is, if i G I q has q(i) = (xi,mi), then the system currently contains a particle of 
mass irti at x j el 2 . To define the initial configuration, we choose a sequence of continuous 
functions {h n : R 2 -> [0, oo), n G N} that must satisfy some conditions that we will shortly 
specify. We set Z = YlnLi fm?hn £ (0, oo), and choose N points in N x R 2 indepedently 
according to a law whose density at (x,n) is equal to h n (x)/Z. Selecting arbitrarily a set of 
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N symbols {ij : j G {1, . . . , N}} from J, we define the initial configuration q Q by insisting 
that qo(ij) is equal to the j-th of the randomly chosen members of NxR 2 . 

Each particle moves according to an independent Brownian motion whose diffusion rate 
2d(m) depends on its mass m G N. As we will explain later, we require some conditions on 
the choice of the function d : N — > R, although the restriction imposed by these conditions 
is far from prohibiting the physically reasonable choice where d is decreasing. Any pair of 
particles that approach to within a certain range of interaction are liable to coagulate, at 
which time, they disappear from the system, to be replaced by a particle whose mass is equal 
to the sum of the colliding particles, and whose location is at some point nearby the place 
where the collision took place. This range of interaction is taken to be equal to a parameter 
e, whose dependence on the total particle number N must be stipulated. We make the 
choice N = \loge\Z . This will ensure that a particle randomly chosen from those initally 
present experiences an expected number of collisions in a given unit of time that remains 
bounded away from zero and oo as N is taken to be high. The effects of motion and reaction 
determine the macroscopic evolution of the system to comparable extents in this scaling. 

We now describe the mathematical details of these dynamics. Let F : {M 2 x N} 1 — > K. 
denote a smooth function, whose domain is given the product topology. The dynamics is 
such that the action on F of the model's infinitesimal generator L is given by 

(hF)(q)=A F(q)+A c F(q), 

where the diffusion and collision operators are given by 

A F(g) = ^d(m i )A Xi F (1.4) 
iei q 

and 

A c F(q) = ^r 2 |logefV(^)«K mj ) (1.5) 

-J^—F^q) + -J^—F^ q) - F(q) . 
rrii + rrij ,J m, + rrij ,J 

Note that: 

• the function V : 1R 2 — > [0, oo) is assumed to be Holder continuous of compact support, 
and with J R2 V{x)dx = 1. 

• we denote by S^q that configuration formed from q by removing the indices % and j 
from I q , and adding a new index from I to which S} fl assigns the value (xj, mi + rrij) . 
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The configuration Sf^q is defined in the same way, except that it assigns the value 
(xj, rrii + rrij) to the new index. The specifics of the collision event then are that the 
new particle appears in one of the locations of the two particles being removed, with 
the choice being made randomly with weights proportional to the mass of the two 
colliding particles. 

We will denote by the measure on functions from t G [0, oo) to the configurations 
determined by the process at time t. Its expectation will be denoted Ejv- 

The form of the collision term in (|1.5J) differs from that used in the case of higher dimen- 
sions, in that the factor of | loge| is absent in the latter case. To explain why we make 
this change, we firstly recall the reason for the form of the collision operator in the case 
when d > 3. Suppose that, for some such choice of the dimension, two particles (xi, rrii) and 
(xj, rrij) have, at some time t , just become liable to interact, in the sense that the difference 
Xi — Xj has become of order e . This state of affairs is liable to persist for a time of order e 2 , 
but not much longer: for d > 3 and C a large constant, the Brownian displacement Xi — Xj 
would return a distance of e from the origin with only a small probability after a time of Ce 2 
after the moment to- This means that, by choosing a form of collision dynamics in which 
the factor of | loge| _1 is absent from (jl.5j) . we ensure that the integral 

f T x x 

It — a(mi,mj)e~ 2 V ( — -Jdt, 

Jt=t V e / 

reaches its eventual value after a time of order e 2 . Other particles are unlikely to interfere 
with this pair in such a short period of time, and, as such, we may neglect their influence. 
The probability of collision between the pair before time T is equal to 1 — exp { — It] ■ 
Thus, for d > 3, our choice of dynamics is such that, among all the pairs of particles that at 
some moment lie within e of each other, the fraction that eventually coagulate is bounded 
in N away from and 1 , with this fraction being close to or 1 depending on whether the 
relevant constant a(mi,mj) is high or low. 

Turning to the planar case, note firstly that, in order that the probability of pair collision 
may remain equal to 1 — exp { — It} , we alter the definition of It by introducing a factor 
of | loge| . The two-dimensional case differs, because a planar Brownian motion returns 
almost surely to any open set at indefinitely later times. As such, the difference x« — Xj will 
endlessly re-enter the e-ball centred at the origin, ensuring that It — > oo as T — > oo. In 
a system of two particles, their coagulation is inevitable. In the system that we consider, 
where the regime of constant mean free path has been selected, a pair of particles at e 
distance may find that their ongoing efforts to coagulate, as measured by the increase of It, 
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are interrupted by the arrival of a third particle, the probability of appearance of such an 
intruder becoming appreciable at a small time, independent of e, after that at which the 
pair in question first came close to each other. The factor of | log el 1 that appears in (jl.5|) 
ensures that, during this short fixed time, the probability of coagulation between the pair 
is of unit order, with the constant a(jrii,mj) determining whether this probability is high 
or low. To write a statement analogous to that for the higher dimensional case: among the 
set of pairs of particles that are at some moment at a distance of order e, the fraction that 
combine with each other, rather than with some other particles, is bounded in N away from 
and 1 , with the value of the constant a determining whether this fraction is high or low, 
similarly to the earlier case. 

Our main result is conveniently expressed in terms of the empirical measures on the 
locations of particles of a given mass. For each neN and t G [0, oo) , we write g n {dx, t) for 
the measure on M. 2 given by 

g n (dx,t) = llog^ 1 2j 6 Xi (t)(dx)l(mi(t) = n). 



id) 



We also require a mild hypothesis on the diffusion coefficients d : N — > (0, oo) (see the 
first remark after Theorem 1.1 below). Namely, we suppose that there exists a function 
7 : N 2 — > (0, oo) such that a < 7, with 7 satisfying 

7127(711, n a + n 3 ) maxjl, rf (^2 + w 3) } < (n 2 + n 3 )-f(n l ,n 2 ). (1.6) 



d(n 2 ) 

The initial random configuration of N particles is formed by scattering particles of nu- 
merous masses independently in M? according to densities that are prescribed for each mass. 
These densities will be chosen as continuous functions {h n : M 2 — > [0, 00), n G N}, and 
should satisfy some fairly weak bounds. To be specific, we insist that 

• fee L X (M 2 ) and k G L™ C {R 2 ) where k := Y,n=i nh n and 



k(x) = // I log \xi — x 2 — y\ \k(xi)k(x2)dxidx2- 

J J \x\ — X'2 — J/|<1 

• For every m, Yln=i d(n) 2 ^ 3 'y(n, m)h n G L^ C {R 2 ) for h n (x) = J h n (y)\x — y\~ 4 ^ 3 dy. 

• For every m, Yl™=i d(n) 3 ^'j(n, m)h n G L^ C (R 2 ) for h n (x) = f h n (y)\x — y\~ 3 ^ 2 dy . 

We then set Z = Y^n=i Ik 2 ^ n e 00 ) anc ^ choose N points in N x I 2 independently 
according to a law whose density at (x,n) is equal to h n (x)/Z. Selecting arbitrarily a set of 
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N symbols {ij : j G {1, . . . , N}} from J, we define the initial configuration go by insisting 
that qo(ij) is equal to the j-th of the randomly chosen members of NxR 2 . 
Remark It is not hard to show that our assumptions on the initial data {h n } are satisfied if 
k is bounded, has a bounded support, d(-) is bounded and 7(71, m) < C{m)n for a function 
C(-). Indeed if k is bounded and has a bounded support, then kL^ c and h n ,h n G L^ c for 
every n. It is worth mentioning that if k belongs to the negative Sobolev Space H~ L = 
W- 1 ' 2 , then keL°°. 

The main theorem is now stated. 

Theorem 1.1 Let J : R 2 x [0, 00) — > K be a bounded and continuous test function. Then, 
for each n6N and t G (0, 00) , 

lim sup Ejv 

where we recall that e is related to N by means of the formula N\ loge| _1 = Z , with the 
constant Z G (0, 00) being given by the expression Z = ^ ngN / R2 h n . In (1.7), {f n : M 2 x 
[0, 00) — > [0, 00), ?i6N} denotes a weak solution to the system of partial differential equations 
(1.1) with the initial data / n (-,0) = h n (-) . The collection of constants (3 : N 2 — > [0, 00) is 
given by 

2n(d(n) + d(m))a(n, m) 

p(n,m) = — -. — -v^ -. (1.8) 

2n\d{n) + d[m)j + a(n, m) 

Remarks 



J {x, n, t) (g n (dx, t) - f n (x, t)dx) 



0. 



;i.7) 



• Conditions for uniqueness among weak solutions of (1.1) are provided in Proposition 
2.6 of Although these conditions are slightly stronger than those we have proved 
the macroscopic equations to satisfy in our case, we anticipate that this uniqueness 
holds. In the lack of uniqueness, the best that can be said is that if Vn denotes the 
probability law of g(dx,t) with respect to the underlying probability distribution Pjv, 
then {Vn} is tight and any limit point V of the sequence \Vn\ is concentrated on 
the space of measures g n (dx,t) = f n (x,t)dt with {f n : n G N} solving (1.1) in weak 
sense, subject to the initial condition /„(-,0) = h n (-). 

• Included in the space of parameter values that satisfy (jl.6|) is the case where the dif- 
fusion rate d is a decreasing function of the mass, and the coagulation propensities a 
satisfy a(n,rn) < Cnm. In fact for a nonincreasing d(-), the condition (jl.fij) is equiv- 
alent to saying that a(n,m) < C(ri)m for a function C(n). Also, if the microscopic 
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coagulation rate a is identically constant, then the condition (jl.6J) is equivalent to 
saying that the function d(n)n~ l l A is nonincreasing. 

• Note that the macroscopic coagulation propensities f3 depend only on the total inte- 
gral of V that is assumed to be 1 for convenience. However when the dimension 
is 3 or more the propensity /3(n,m) does depend on V in a nontrivial way and 
is given as a(n,m) J(l + u)Vdx, where u solves the PDE Am = r(l + u)V with 
t = a(n, m) / (d(n) + d{m)) . 

• Our technique of proof also yields a kinetic limit derivation for the model in which 
particles are assumed to have a range of interaction that is mass-dependent. To give 
an example of such a variant, suppose that each particle of mass m has a radius r(m) , 
where rim) = y/m. We stipulate that particles of mass m and n are liable to react 
when their displacement reaches the order of (r(m)+r(n))e. More precisely, we modify 
the definition (jl.5|) of the collision operator Aq by replacing the appearance of V by 
(r{n) +r(m)) V( ■ / (r(n) +r(m))) , (the factor that multiplies V being introduced so 
that, roughly speaking, the altered collision mechanism respects the spatial-temporal 
scaling of Brownian motion). Theorem 1.1 is still valid for this modified model with the 
same macroscopic coagulation propensities (3 . This is in sharp contrast with the case 
d > 3 for which the mass dependence affects the macroscopic coagulation propensities 
/3. 

In common with the proof for d > 3 , a central element in deriving Theorem 11.11 is 
establishing that, at any given moment after the initial time, the presence of a particle of 
some given mass at some fixed point in space significantly affects the likelihood of a particle 
being at some other point in space only if that other point is at a short distance from the first 
particle. That is, on distances of short order, the presence of a particle makes it less likely to 
find another nearby, because the pair would have been liable to coagulate shortly beforehand. 
However, the distribution of particle at a given time is similar to one in which they were 
scattered independently, except for this short-range repulsion. The following proposition, 
whose form differs from that in the case d > 3 only in its scaling factor, formalises this 
assertion. 

Proposition 1 Set 

Q = \\oge\~ 2 a(mi,mj)V t {xi — Xj)J{xi,rni,t)J{xj,rrij,t), (1.9) 

where J, J : M 2 x N x [0, oo) — > [0, oo) are test functions satisfying the same conditions 
as those stated in Theorem li.il We also assume that J(x, m,t) =0 unless m = Mi and 
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J(x, m,t) = unless m = M 2 . Let rj : M 2 — > [0, 00) denote a smooth function of compact 
support for which L 2 7)(x)dx = 1 . We have that 



[ Q(t)dt 
Jo 

I dt I dcuP(M 1 ,M 2 )J(uj,M 1 ,t)J(u,M 2 ,t) \\oge\- 1 Yl <S "'/('' 



1.10) 



X.j — LU 



5 

j£l q ;mj=M2 



+ Err(e,5), 



where the constants (3 : N 2 — > [0, 00) were defined in hi. 8$ . and where the function Err 
satisfies 

lim lim sup Ew I Err(e, 5) \ = 0. 
sio 40 1 v 1 

Why is this statement a mathematical rendering of the claim discussed before it was made? 
The quantity f Q Q(t)dt can be thought of as the total propensity of particles to combine 
during the interval of time [0, T] . Proposition 1 asserts it may be approximated by a time- 
averaged product of empiricial approximations to the density of particles (of the appropriate 
mass) . That is, particles are arranged independently enough near most of the collision events 
that the rate of these collisions is roughly proportional to that arising in a system in which 
particles are scattered indepedently at random according to densities given by measuring 
the system in question on scale 5 that is much larger than the reaction range e. There 
is, however, a constant of proportion corresponding to the change from microscopic reaction 
propensity a appearing in the definition of Q to its macroscopic counterpart (3. Its presence 
may be explained by the negative short-range correlation between particles discussed before 
the statement of Proposition [TJ 

The analogue of Theorem 11.11 that appears in [2] for the case d > 3 is derived as a 
consequence of Proposition ^ In Section 2 of |2j, a sketch of the proof of Theorem 11.11 mav 
be found. The details of the derivation of Theorem ll.ll from Proposition^ do not differ in the 
two-dimensional case, so that we do not present these arguments again in this paper. Our 
task here is rather to present a detailed derivation of Proposition ^ in the case when d = 2. 
Before reading further, however, the reader may wish to consult Section 2 of We refer the 
reader to [2] also for a discussion of previous work related to the problem. Here, we mention 
only Sznitman 4], in which a model of Brownian spheres that annihilate as soon as they 
touch is studied. The partial differential equation by which the density of particles evolves 
was derived for the kinetic limit, in each dimension d > 2. In this work, the macroscopic 
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annihilation rate is exactly 2ir when the dimension is 2. This is compatible with our main 
results because if d(-) is identically 1/2 and a — > oo, the macroscopic coagulation rate (3 
approches 2tt . Note that our model approximates the hard core model as a gets large. 



2 Establishing the Stosszahlansatz 



For any given pair (n, m) 6 N 2 , w e = m : M 2 — > [0, oo) is a function whose existence is 
ensured by Theorem 13. II that lies in C 2 (lR 2 ) satisfying 



(2.1) 



<,mW = ^ r ( n > m ) J - y\[V e (y)u e niTn (y) +V e (y) 

where r(n,m) = a(n,m)/(d(n) + d{m)). As a consequence, 

(d(n) +d(m))Au £ nm (x) = a(n,m) V e (x)u e n>m (x) +V e (x) 
We are using the notations 



and 



^ e (x) = e-V(- 



V e (x) = e- 2 |loge|-V( 



We present the conditions on the two test functions J, J : M 2 x N x [0, oo) — > R that 
appear in Proposition 1. It suffices to work with functions that take non-zero values for only 
one value in the second argument, such functions measuring the presence of particles of a 
given mass. By a temporary abuse of notation, we write 

J(x,M u t) = J(x,t)l{m = Mi} 
J(x,M 1 ,t) = J(x,t)t{m = M 2 }, 

where on the right- hand- side, J and J denote smooth maps from M 2 x [0, oo) to R of compact 
support. We will suppress the appearance of the t-variable when writing the arguments of 
J and J. 

In seeking to verify Stosszahlansatz, we define 
X z {q) = | loge|~ 2 u li u M 2 {.Xi-Xj + z)J{x i ,M 1 ,t)J{x j ,M 2 ,t)\{m i = M^rrij = M 2 }. 

(2.2) 
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The relevance of the expression (2.2) for our purposes is that the term Q and its variations 
appear as we apply the infinitesimal generaor on the expression X z — Xq. We refer the reader 
to Section 2 of j2] for some heuristic justification of the special form of X z . 

Numerous terms arise when the operators A and A<^ act on the expression X z — X Q 
(recall that the functions of configurations X z , indexed by z e R 2 , were defined in (2.2)). 
We now label these terms. Unless stated otherwise, we will adopt a notation whereby all 
the index labels appearing in sums should be taken to be distinct. This includes the case of 
multiple sums. For example, ^2 k i Gl Yuiei f( x k, x h x i) denotes the sum of the evaluation of 
the function / over all arguments that are triples (xk,xi,Xi) where k, I and i are distinct 
indices in I . Note also that, unless otherwise stated, whenever the symbol u e appears in a 
summand, we mean u e M M . 

Firstly, we label those terms arising from the action of the diffusion operator. To do so, 
note that, for a time-dependent functional F of the configuration space, this action is given 
by 



(^ t +^)F = ^F + J2d(m t )A Xt F 

iEla 



Thus, we label as follows: 



with 

H 12 
#13 
H\4 



0_ 

Of 



+ A {X z - X )(q(t)) = H n + H 12 + H 13 + H u + H 2 + H 3 + tf 4 , 



loge 



-2 



V £ (%i — Xj + zj — V e (xi — Xj) J(xi,m i ,t)J(xj,mj,t) 



a(m,i, rrij) 

— | log e | 2 a(mi,mj)V e (xi — Xj)u e (xi — Xj)j(xi,mi,t)J(xj,mj,t) 
| log e | 2 a(mi,mj)V e (xi — Xj + z)u e (xi — Xj + z)J(xi,mi,t)J(xj,mj,t), 
| log e| 2 jw e (xi — Xj + zj — u e ( 



X i X j ) 



[j t (x i ,mi,t)j(xj,m j ,t) + J(xi,m i , t) J t (xj , rrij , t) 
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along with 



H 2 = 2 1 log e | 2 d(mj)J(xj,mj,t) 
i,jei q 

^U x (Xi Xj -\- z) U x (^Xi %j) ' Jxip^ii TTbij t) j 

H 3 = — 2| log e| 2 d(m j )J(x i ,m h t) 

^u e x (xi - xj + z) -u e x (xi - xj) ■ J x (xj ,mj,t), 



and 
H 4 



loge| " — Xj + z) — u £ (xi — Xj) 

i,j€l 



^d(mi)A x J(xi, rrii, t)J(xj,rrij, t) + d(rrij)J(xi, mi, t)A x J(xj,m,j,t) 

where f x denotes the gradient of /, and ■ the scalar product. As for those terms arising 
from the action of the collision operator, 

A C (X Z - X )(q) = G z (l) + G z (2) - G (l) - G (2), 

where G z (l) is set equal to 

^2 a(m k ,mi)V e (x k - xi)\\oge\ 

k,lGl q ielq 



m k + mi 



+ 



m k + mi 



u\x k - Xi + z)J(x k ,m k + m l ,t)J(x i ,m i ,t) 

+ u e (xi -x k + z)J(xi,mi,t)J(x k ,m k + m h t) 
u e (xi -Xi + z)J(xi,m k + mi, t)J(xi, m u t) 

+ u e (xi -xi + z)J(xi,mi,t)J(xi,m k + m h t) 
u e (x k - Xi + z)J(x k ,m k ,t)J(xi,mi,t) 

+ u e (xi - x k + z)J(xi,mi,t)J(x k ,m k ,t) 
u e (xi - Xi + z)J(xi, m x , t)J(xi, mi, t) 

+ u e (xi -xi + z)J(xi,mi,t)J(xi,mi,t) 
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and where 



G*(2) 



loge| 2 ^2 at(m k ,mi)V e (x k - x\) 
k,iei q 

u e (x k - xi + z)J(x k ,m k ,t)J(xi,rn l ,t). 



(2.3) 



The terms in G 2 (l) arise from the changes in the functional X z when a collision occurs due 
to the influence of the appearance and disppearance of particles on other particles that are 
not directly involved. Those in G z (2) are due to the absence after collision of the summand 
in X z indexed by the colliding particles. 
Note that 

H 12 + G (2) =0. (2.4) 
The process { (X - X ) (t) : t > 0} satisfies 



(- + A ) (X z -X )(t)dt 



(2.5) 



+ / A c {X z -X )(t)dt + M(T), 



with {M(t) : t > 0} being a martingale. By using the labels for the various terms that we 
just introduced, we find from (|2.5jl by use of (|2.4jl that 



f 

H n (t)dt+ / H n (t)dt 
Jo 

< \X z -X \(q(T)) + \X z -X Q \(q(0)) 



+ 



r^~i rj~\ r^~i 

[ \H u \(t)dt+ [ \H 2 \(t)dt+ [ \H 3 \(t)dt+ [ \H 4 \(t)dt 
Jo Jo Jo Jo 



+ / \G Z (1)-G (l)\{t)dt+ \G z (2)\(t)dt + \M{T) 



(2.6) 



Since J is of compact support, we have that X z (q(T)) = for T sufficiently large. We aim 



12 



to prove the following estimates: for each T > 0. 

f -T 



E N \H 1A \(t)dt < Clz^llog 





fT 



[ E N \H 2 \(t)dt<C\z\ 1/9 \log\z\\, 
Jo 

[ E N \H 3 \(t)dt < C\z\ 1/9 \\og\z\\, 
Jo 

I E N \H A \(t)dt<C\z\ l/2 \\og\z\\, 
Jo 

[ E N \G Z (1) - G (l)\(t)dt < C\z\ l ' 2 \ \og\z\\\, 
Jo 

[ E N \G z (2)\(t)dt < C\ log |^|| | logel" 1 , 
Jo 

E N \X z -X {Q)\<C\z\. (2.7) 

Later, we apply the limit \z\ — > after sending e to 0. We will also show that, for each 
T G (0, oo) , 

Ejv[m(T) 2 ] < C\ logep 1 . (2.8) 
2.1 Lemmas bounding collision propensity 

In this subsection, we discuss three lemmas that in essence serve as the backbone of the proof 
of the various inequalities that appear in (2.7). These lemmas allow us to reduce the proof 
to a calculation involving the initial configuarions for which the independence of particles 
and our assumptions on the initial densities can be used. In fact the proof of Lemmas 2.1 
and 2.3 is very similar to the corresponding Lemmas 3.1 and 3.3 of For this reason, their 
proofs are omitted. It is Lemma 2.2 that is somewhat different from what we have in [2] as 
Lemma 3.2 and we provide a detailed proof for it. In fact this difference explains to some 
extent a major technical difficulty that is two dimensional and is not encountered when the 
dimension is 3 or more. To explain this further, let us observe that if the dimension d is 3 or 
more and J is a nonnegative function, then we can find a solution to the Poisson equation 
—AH = J that satisfies H > . Indeed the solution H is defined by 

c (gZ) / \x - y\ 2 ~ d J(y)dy, 
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where c (d) = (d(d — 2)u(d)) 1 with uj(d) denoting the volume of the unit ball in R d . This 
is no longer true in dimension 2 because the solution is given by 



J \og\x-y\ J(y)dy. 



This causes some difficulty in treating various terms that appear in (2.7). To get around 
this, let us define 

H(x) = / log \x - y\ J(y)dy. (2.9) 

Zn J\x-y\<l 



We now have that — AH = J — J where 



J( x ) = -Ll J( x -z)dS{z), (2.10) 
' '|*|=i 



2tt 



where dS denotes the 1-Lebesgue measure on the unit circle S 1 . The point is that by using 
Lemma 2.2, we reduce bounding an expression involving J to an expression involving H at 
time t — 0, and a similar expression involving J. Since the funcion J is an average of J, 
we have an easier task to bound the expression involving J . In the case of the terms H 2 
and H 3 , we need to apply this process three times so that the final J has a simple pointwise 
bound. Our three lemmas are: 

Lemma 2.1 For any T G [0, oo) 7 

| log e | Ejy / dt a{mi,mj)V € (xi — Xj) < Z. 

Lemma 2.2 Let J : R 2 — > [0, oo) be continuous, and let H : IR 2 — > [0, oo) be given by (2.9) 
We also define J : M 2 — > R according to (2.10). Then we have the following inequality, 

2 f T 

| loge | / J(xj — Xj)m,imj(d(mi) + d{rrij))dt 



+ 1 log e | 2 Ejy / y ] V £ (xj — Xj)a(mi,mj)H(xi — Xj)mjmjdt 
< | log e | 2 Ejy H{xi — Xj)rriimj 

+ | loge | K N / J(x,i — x j)ra i mj{d{m, i ) + d(m,j))dt. 



f E K(= 
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Lemma 2.3 Assume that the function 7 : N 2 — > (0, 00) satisfies 

i \ f r^( n 2 + ^3)1 2 i / \ 1 \ 

n 2l ( ni ,n 2 + n 3 ) max(l,[ — — J j < (« 2 + n 3 )j(ni, n 2 ), (2.11) 

There exists a collection of constants C : N 2 — > (0, 00) , such that, for any smooth function 
J : IR 4 — > [0, 00) , anc? an?/ gii> en n 1; n3 G M, 

Eat/ <it 2j 7( m i,^j)K(^ - Xj) J(a; i ,a; fc )l{m i = n 1? m fc = n 3 } 

where, also given e > and n 2 G N, the function A% 11 n2 n3 : IR 6 — > [0, 00) is defined by 

(d(ni)A Xl + d(n 2 )A X2 + d(n 3 )A X3 )A e nun2n3 (x 1 ,x 2 , x 3 ) 

= -j(n 1 ,n 2 )V e (x 1 - x 2 )J(x 1 ,x 3 ). (2.13) 

It is worth mentioning that the function n2 n3 (xi, x 2 , x 3 ) of Lemma 2.3 is given by 

Co(6) / / / j(n 1 ,n 2 )J(z,y')V e (z-y)dzdydy', 
Jn? 7r2 V rf(ni) d(n 2 ) d(n 3 ) / 

where co(rf) = (d(d — 2)u(d))~ 1 with c<j(d) denoting the volume of the unit ball in IR d . Note 
that for Lemma 2.3 we are dealing with a solution to a Laplace type equation in M 6 as 
opposed to Lemma 2.2 for which the pecularity of the Laplace equation in M 2 played a role. 
This is why the proof of |2| in the case of Lemma 2.3 can be repeated line by line. 
Proof of Lemma 12.21 Set 

X q = J loge| 2 H(xi — Xj)mimj. 

Recall the mechanism of the dynamics at collision: the location of the newly created particle 
is one of the two locations of the colliding particles, with weights proportional to the masses 
of the incident particles. We see that when Ac acts on X q , all those terms indexed by pairs 
of particles one of which is not involved in the collision cancel. Thus, 

1 1 —2 X — \ 

A C X = — J loge| y, V<i\ x i ~ Xj)a{mi,mj)mirajH{xi — Xj). (2-14) 
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By AH = -J + J, 

I 1—2 X — v 

A X = loge y j AH[Xj — Xj)mimj\d{mj) + dyrrij)) (2-15) 

I I — 2 X — * 

= — | log e | y j J{x; L — x j)m i mj{d{m, l ) + d(m,j)) 



-2 



I I— 2 X — ^ ~ 

+ | loge| " y^ J{Xi — Xj)mimj(d(mi) + d(rrij)) 

= : AlX + A 2 X (2.16) 

From the non-positivity of A^X , the non-positivity of AqX , apparent from (|2.14|) . and the 
non- negativity of X , follows 

-E N f AlX(t)dt-E N [ A C X (t)dt <E N X(0) +E N [ A 2 X(t)dt. (2.17) 
Jo Jo Jo 

□ 



2.2 Bounds on functionals of u. 



We will verify the assertions presented in (|2.7jh The following lemma provides the bounds 
on the behaviour of the functions {u e nm : M 2 — > [0, oo) : (n,m) G N} and other functions 
that will be used in this section. We choose the constant R so that V(x) = whenever 
|x| > R . Recall that k = ^2 n nh n . 

Lemma 2.4 There exists a collection of constants C : N 2 — > (0, oo) for which the following 
bounds hold. 



for x G R 2 satisfying \x\ < 2R e , \u e n ^ m (x)\ < C niTn |log£| ; and for all x G 
|< >m (z)| < Cn,m| log | a: 1 1 . 

for x G R 2 , |V< im (x)| < C„, m min|^,i}. 

for x <ER 2 , 



< Cr, m I z I min 



1 1 

\x\ ' e 



(2.18) 
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for i£l 2 satisfying \x\ > max {2\z\ +R e, 2R e}, 



< 





2 






2 



(2.19) 



Ze£ = i/„ im : I 2 x I 2 ^ [0, oo) be given by 



H(x; z) 



-1 

2tt 



\x-y\<l 



\og\x - y\u e njTn (y + z)l{\y\ < p}dy. 



Then, 



sup / H(x x - x 2 ; z)k(xi)k(x 2 )dx 1 dx2 < C n)m p 2 \ logp|. 

\z\<lJ 



Let H( = H njjn ) : R 2 x R 2 -> [0, oo) 6e owen 6y 



H ( x ; z )=2^[ log\x-y\ Vu e n>m (y + z) l{\y\ < p}dy. 



\x-y\<l 



Then, for every z with \z\ < 1, 

J H(xi - x 2 ; z)k(x 1 )k(x 2 )dx 1 dx 2 < C n>m (p + \z\). 

for (x, z) G IR 2 x M. 2 , let L(x; z) be given by 



(2.20) 



log|x-y| [ | log j 1 — \y + z\\\ + 1] l{l-p< + < l + p}dy. 

Zn -/|x-j/|<1 



T/ien we nave tne bound 



L(x;z) < Cp(logp) 2 . 



(2.21) 



/or any positive integers n and m and a nonnegative smooth function J of compact 
support, there exists a constant C n<m (J) such that, for any given z6K 2 , the function 
■> [0, oo) defined by 

(d(ni)A xl +rf(n 2 )A :E2 + d(n 3 ) A a . 3 )A^ lin2in3 (xi, x 2 , x 3 ) 

* >m (ici - x 3 + ^^(xi - x 2 )l{|xi - x 3 \ < p}J(x 3 ) 



satisfies 



52 / A t nin2na (x 1 ,x 2 ,x 3 )h ni {x 1 )h n2 (x 2 )h n3 (x 3 )dx l dx 2 dx 3 

< C n , m (7)| loge| _1 (p + |z|) log (p + |z|). (2.22) 



m,n 2 ,n 3 
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Proof Throughout the proof, we write u £ for the function u e n m and r for the constant 
a(n,m)/(d(n) + d(m)) . The dependence of the constants on n and m arises from that of 
r , and is also omitted. The first part of the Lemma is a straightforward consequence of our 
results in Section 3. As a consequence of Theorems 3.1-3.2 and Lemma 3.1 we know that 
there exists a constant C\ G (0, 1) such that for small e and y satisfying \y\ < 2R , 

loge < u £ (ey) < c\ logs, 

or equivalently, 

< u e (ej/)|loge| _1 + 1 < 1 -ci. (2.23) 
From (2.1), we learn that for x satisfying \x\ > 2eR , 

A £ \og(\x\/2) < u £ (x) <\og(2\x\)A £ , 

where 



A £ = ^r| (y e {y)u* n>m {y) + V%y))dy. 



From this and (2.22) we learn that there are two positive constants k\ and k 2 such that if 
\x\ > 2eR , then 

fcilog(|a;|/2)) <u £ (x) < fc 2 log(2|a;|). 
To prove the second part of the lemma, recall firstly that 

u e (x) = ^ J 2 lo £ \ x - v\ (u e (y) v e(y) + ve (y)) d y- 

As a result, 

Z7T Jm. 2 \x - y\ v ' 
If |x| > 2R e, then \x — y\ > \x\/2, and 



x-y 



\x - y\ 



1 2 

r < — 



\x - y\ 



implying that 

fl \ i 2A e 3 
V« £ x < -p-f < — , 
1 1 \x\ \x\ 

for small e. If \x\ < 2R e, then we use (2.22) to deduce 



r i rittQt 
\S7u e (x)\ < Cl / , TV £ (y)dy < c 2 e~ 2 / dr < c 3 ^ ] 

J\x-y\<3Roe \ x ~ 2/1 ./0 



18 



Thus, 



Vw e (:r) < C min 



1 1 

e \x\ 



as claimed in the second part of the lemma. 

To prove the third part of the lemma, note that 



u 



[x + z) — u e (x) 



< — 

1 

< — 



log 
log 



\x — y + z\ 



\x - y\ 
\x — y + z\ 



\x - y\ 



u e (y)V e (y) + V £ (y) )dy 
V € (y)dy, 



(2.25) 



the latter inequality by means of (2.22). From this and the elementary inequalities 



log 



\x — y + z\ 



we deduce that 



x -y\ 
u e (x + z) — u e {x) 



1 < 



\x — y + z\ 
\x - y\ 



< 



< 



2tt 



1 

27 



\x - y\ 

v € (y) 

\x - y\ 



(2.26) 



dy 



We now use this and argue as in the proof of the second part of the lemma to deduce the 
third part of the lemma. 

In seeking to prove the fourth part of the lemma, note that 



x + z — y 




x-y 




x-y 


2 (x + z - y) - 


X + z ■ 


-y 


2 {x-y) 


x + z — y 


2 


x-y 


2 


x + z — y 


2 


\x-y 


2 



Note that, for any o£l 2 , 

\a\ 2 (a + z) — \a + 



z\ 2 a 



< \a 

< Ci\z\\a 



a + z\ — \a\ 

2 



+ \z\\a\ 



(2.27) 



so long as \z\ < \a\. Note that since by our assumption |x| > 2\z\ + RqE , we have that 
\x — y\ > 2\z\. We may apply (|2.27|) with the choice a = x — y to the formula (|2.24|) . hereby 
obtaining 



Vu e {x + z) - Vu e {x) <C\z\ [ 

Jr 2 



v e (y) 

\x + z — y\ 



,2 d yi 
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where we used (2.22). From the inequality \x\ > max{2|z| + R e, 2it! e} , we deduce that 
\x + z — y\ > \x — y |/2 and \x — y\ > \x\/2. We conclude that 



Vu £ (x + z) - Vm £ (i) 



as required. 

To prove the fifth part of the lemma, note that 



x 



H(x;z)\ <C |log|a;-y||log|y + z||l{|y| < p}dy , 

J\x-v\<l 



'\x-y\<l 

by the first part of the lemma. Hence, 



J \H(xi — x-i\ z)\k{xi)k{x2)dx\dx2 

log \xi — x 2 — y\\k(x l )k(x 2 )dx 1 dx 2 dy 



< 


c 


j 






l\y\<p 


< 


c 


! 






)\y\<p 








< 


c 


r i 






lo 



' \xi-x2-y\<l 



\og\y\\dy 



where in the second inequality, we used our first assumption on the initial data and the fact 
that if p + \z\ < 1, then the expression f^ <p \ \ log \y + z\\dy\ is maximized as a function of 
z G 1R 2 when z — 0. This establishes the third part of the lemma. 
To prove the sixth part of the lemma, note that, by the first part, 

1 f 1 

H(x;z)<— \\og\x-y\\- — ■ — -l{ \y\ < p}dy. 
2tt J\ x -y\<i \y + z\ 



It follows that 



J \H(xi — x 2 ; z)\k(xi)k(x 2 )dxidx 2 

< ci / | — ^ — f / \log \x! - x 2 - y\\k(x 1 )k(x 2 )dx 1 dx 2 dy 

J\y\<p \V + A 7|xi-x 2 -S/|<l 

< c 2 / | — ■ — -dy <c 2 - — ■ — -dy 

J\y\<p W + Z\ J\y+z\< P +\z\ \y + z\ 

< c 3 (p+\z\). (2.28) 
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We have deduced (l2~2nj) . 

As for the seventh part of the lemma, first observe that L(x\ z) = L(x + z; 0) . Hence we 
only need to verify (|2.21|) when z = 0. In this case we divide the domain of integration into 
the sets \x — y\ < |1 — \y\\ and |1 — \y\ \ < \x — y\ < 1. Hence, 



L(.;0) < - 



\*-y\<p 



ii 1 2 

log \x — y\ + log \x — y\\ dy 



1 

+ /I 

Z7F ^ 1— Is/I \<\x-y\<l 



log |1 - b||| + | log |1 - bHl 2 G (1 — p, 1 + p))dy 



< I r 
'o 

< C'p(logp)^ 



i+p 



[| logr| + | log r| 2 ] dr + / | log |1 — r\ | + | log |1 — r 



i-p 



rdr 



establishing f[2~2"T|). 

As for the eighth part of the lemma, let us write J(a) for l{|a| < p} and define the 
quantity / according to 



/ = cq(6) 7(^1,^2) \dx\dx 2 dx 



\x\ - z'\ 2 \x 2 - y\ 2 \x 3 - y'\ 2 ^ ~ 2 



(ifrii) (if^?) cifr^) 

ni,ri2,n3 " " v/ v/ v/ 

/i ni (xi)/i„ 2 (x2)/i„ 3 (x3)w £ (z / - y'+ ^KfV - y)J{z - y')J(y')dz'dydy'. 

We write 

J = c (6) / u%3> -tf + z)V& -y)J{J -rtlWGy^yWdydtf, 

J MP 

where G(z',y,y') is given by 



ni, 7i2, 7-137(711 ,712) 



H 77 — ; 1 — — ; — I h ni (x 1 )h n2 (x 2 )h n3 (x 3 )dx 1 dx 2 dx 3 . 



1-4/3 



d(ni) d{ri2) d{n 3 ) 

Using the elementary inequality abc < (a 2 + b 2 + c 2 ) 3 / 2 we deduce that G(z', y, y') is at most 

^ 7(ni, n 2 ){d(n 1 )d(n 2 )d(n 3 )f l ' i J \x x - z'\~ 4/3 \x 2 - y\~ 4/3 \x 3 - y'\ 

h ni (xi)h n2 (x 2 )h n3 (x 3 )dxidx 2 dx 3 . 
From our assumptions on h n we deduce that G G Lf^ c . Hence, 

I<C [ u £ (z' -y' + z)V e (z' - y)J(z' - y')J(y')dz'dydy' . 
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Note that, for fixed z' G 



Jm^ \ e J 



Thus, 



/ < Close 



< Close 



< 




u e {z! -y' + z)J(z' - y')J(y')dz'dy' 
u e {z — y' + z)J(z' — y')dz'dy' 
Cjlogep 1 / dy' [ \\og\z' -y' + z\\j(z' - y')dz' (2.29) 
log \a\\da < C| loge| (p + \z\) 2 \ log (p + \z\)\, 



< Close 



\a\<p+\z\ 



where K C R 2 denotes a compact set containing the support J, and where we made use of 
the first part of the lemma in the third inequality. This is the bound stated in ()2.22|) . □ 

2.3 Estimating the terms 
2.3.1 The case of H u and H 4 

The estimate of E, N f Q T \H 4 (t)\dt is derived in an identical fashion to that of E, N f Q T \Hi 4 (t)\dt. 
Note that 



E 



N 



H 14 (t)dt 



2 f T I 

< C|loge| Ejv / dt \u € (xi — Xj + z) — u € (xi — Xj) 
1 {mi = M^rrij = M 2 }, 



where the constant C depends on the L°° bounds satisifed by J, J and their time derivatives. 
Hence 



E 



N 



H u (t)dt 



<K, + K 2 , 



where K\ is given by 



|-2 T 



C| log e| Ejy / dt 



i,jeI q :\xi-Xj\>p 

l{ mi = M^linij = M 2 } 
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and K 2 is given by 



loge E 



dt u e (xi — Xj + z) — u e (xi — Xj) 

i,jeI q :\xi—Xj\<p 

l{mi = M^lirrij = M 2 }. 



Firstly, we treat K\. Note that 



C\z\ | loge | 
P 



-E 



N 



dt H m i = M^lirrij = M 2 } < 

i,j£l q :\xi-Xj\>p 



C\z\Z 2 
P 



where the first inequality follows from the third part of Lemma 12.41 and the final one from 
the initial number of particles N equals Z\ loge \ . 
We now treat the term K 2 . By writing, 



-2 T 



Ko < C loge E 



N 







U e (x,i — Xj + z) I + \u e (Xi — Xj) 



dt J2 

i,j€lq.\Xi-Xj\<p 

l{ mi = M x }\{mj = M 2 }, 

we obtain an expression on the right-hand-side which may be bounded by applying Lemma 
2.21 As a result we can write K 2 < K 2 \ + -^22 where K 2 \ and K 22 represent the first and 
the second term on the right-hand-side in Lemma 12.21 For K 2 \ , the relevant estimate is 
provided by the fifth part of Lemma l2.4[ with a bound of Cp 2 |logp|. To bound the term 
K 22 , note that, with the function J in Lemma E21 chosen to be J(x) = u e (x + z)1(\x\ < p) , 
we have that 



\J(x)\ 



u 



bl=i 



(x + y + z)\{\x + y\ < p\S(dy) 



< C \\og\x + y + z\\l{\x + y\< p}S(dy) 
J\y\=l 



'\u\ 

Let us assume that \z\ < p ans set a = x + z. We then have 



\J(x)\ < C 



log|a + y||l{|a + y| < 2p}S(dy) 



/•C2p 

< 2C |log(cir)|dr < c 3 p|logp|, 
Jo 



where for the second inequality we used that fact that the conditions \a + y\ < p, \y\ = 1 
mean that the point y belongs to an arc on the unit circle of center a = a/\a\ and length 
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2c 2 p, and if the lenght of the arc ay is r, then \a + y\ > c±r for positive constants c\ and 
c 2 - 

In this way, we find that J is uniformly bounded by Cp\ logp| and this in turn implies 
that the term K 22 is bounded above by 

C\ loge|~ 2 p| logp|E,v J dt\{{i,j) e l\ : m* = M^m, = M 2 }| < CZ 2 p| logp|, 

the latter inequality following from N < Z | log e | . We find that 

X 2 <Cp|logp|. 

Hence, 



E 



#i 3 (*) 



dt < K x + X 2 < + Cp\ logp|. 
P 



Setting p = | z I 1 / 2 , we find that 



E 



A? 



Jo 



dt < Cj^| 2 1 log |^| I . 



2.3.2 The cases of H 2 and 



The estimate of E^r J \H 3 (t)\dt is derived in an identical fashion to that of E^r f \H 2 (t)\dt. 
Picking pGl that satisfies p > max {2\z\+R 6,2R e}, we write 



/ E N \H 2 (t)\dt < Ri + R 2 , 
Jo 



where 



Ri = 



loge|~ 2 Ejv / ^ d{m i )\j{x j ,m j ,t)\ 

° i,j&I q :\x t -Xj\>p 



It, 



and 



#2 



lo K e E 



N 



I ^2 d(mi)\j(c 

i,j&Iq-\Xi-Xj\<p 

u e x (xi - Xj + z) - u e x (x,i - xj) | J x (xi, m h t) \dt. 



(2.30) 
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Firstly, we examine the sum R\ . Recalling that we consider test functions J and J respec- 
tively supported on particles of mass Mi and M 2 , 



1-2 \Z 



Rx < C\ loge| ^d{M x ) e I 2 q : m, = M^m, = M 2 } 

where the lower bound on p allowed us to apply the fourth part of Lemma 12.41 Thus, 



Ri < C^. 



(2.31) 



Secondly, we bound the sum R 2 . Note that 



Ro < lose II J T 



d{Mx) E N J dt 

l{m i = M^rrij = M 2 } 

2 f 1 

< C | log e | Ejv / dt miimj(d(mi) + d(rrij)) 

l{\xi - xj\ < p|, (2.32) 



U^.(^X{ Xj ~\~ Z^j ~\~ U^,(^Xi 



where || ■ || denotes the L°° norm and the constant C depends on the test functions J and 
J. The expression 1)2.32)1 is written in a form to which Lemma 12.21 mav be applied. Doing 
so yields 

R 2 < C|loge| 2 Ejy mivrij {h{xi — xf, 0) + H(xj — %j \ z] 

+ C\ logep 2 E N dtJ2 

J ^ x 2 x j , z ^ I ( x ^ 3j j ^ ^ 



: -R 2 i + R. 
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where the function H appears in the sixth part of Lemma l2~4"l and where J : M 2 xM 2 
in this case is given by 



(2.33) 
[0,oo) 



J{x;z) = — 

ZlT 



\u e x (x + y + z)\t{\x + y\ < p}dS(y). 
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From the sixth part of Lemma 12.41 and our assumptions on the initial data, we deduce 
that 

Rai<C(p+\z\). (2.34) 
It follows from the second part of Lemma f2. 41 that the function J satisfies the bound 

\J(x;z)\<C [ -\{\x + y\< P }dS{y). 

J\ y \=i \x + y + z\ 

By our assumption, we certainly have \z\ < p/2. Hence, 

\J(x;z)\<C [ ^L^l{\a + y\<2p}dS(y), 
J\ y \=i \a + y\ 

for a = x + z. Note that there exists a positive constant c\ such that the conditions 
|« + y\ < 2p, \y\ = 1 mean that t/GT, where T is an arc of the unit circle with the center 
a = — a/\a\. It is not hard to show that there exist positive constants C\ and c 2 such that 

ci(£+ |l - |o||) < \a + y\ < c 2 (£ + |l - |a||), 

where I denotes the length of the arc from a = —a/\a\ to y on the unit circle. From this 
we deduce 

\J{x;z)\ < 1{1 - p < \a\ < 1 +P}J 1 (ciCs 1 ]! - \a\\ < £ < 2qV) jdi 
< c 3 l{l-p< \x + z\ < l + p] ( | log 1 1 — + +log(c 2 q 2 )) 

Using this bound on J , and then applying Lemma 12.21 we learn that the term i? 22 is 
bounded above by i? 2 2i + -R222 , where 



-R221 = C\ loge| 2 E N {L(xi - Xj\ 0) + L(xi - Xj; z) j 



mirrij , 



and, 



R222 = 2C| log e| 2 E N J (^L(xi - Xj) + L(xi - xj + z)^mimj(d(mi) + d(mj))dt, 

where the function L : M? — > [0, 00) appears in the seventh part of Lemma 12.41 and where 
the function L : M 2 — > [0, 00) is given by 

L(x) = — [ l{l - p < \x + y\ < l + p} \\ log |1 - \x + y\\\ + ll dS(y) 
2n J\y\=i 
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By the seventh part of Lemma 12.41 and our assumption on the initial total density k = 
J2 n n hn G L 1 we obtain 

R 221 < Cp(\ogp) 2 . (2.35) 
We decompose R222 = -R2221 + -R2222 where 

R 2 22r = 2| loge| 2 E N J (^L r (xi - Xj) + L r (xi - Xj + z)^mimj(d(mi) + d(mj))dt, 

for r = 1,2, where L\(x) = L(x)1(\x\ > p 1 / 4 ) and L2(x) = L(x)H(\x\ < p 1 ^ 4 ). Let us now 
analyse the behaviour of the function L. First observe that since the Lebsegue measure 
on the circle is rotationally invariant, we have that the function L is radially symmetric. 
Because of this, let us assume that x = (a, 0) and y = (cos9, sin 9) . Note that 

\x + y\ 2 = 1 + a 2 + 2a cos 9. 

As a result, the condition + G (1 — 2p, l + 2p) for p < 1 implies that |a 2 + 2acos$| < 8p. 
Let us first examine the case \a\ = \x\ > p 1//4 . In this case we have that |a + 2cos0| < 8p 3//4 . 
This condition is not satisfied unless \a\ < 2. In that case, choose 9q G [0,7r] such that 
a + 2cos#o = 0. Note that for smallp, the set {9 : \a + 2cos0| < 8p 3 / 4 } is a union of two 
disjoint ^-intervals about the points 9q and 9\ = 2n — 9 Q . We now argue that there exists a 
positive constant C\ such that the length of these intervals is bounded above by Cip 1 / 2 . To 
see this, first observe that the condition \a\ > p 1 / 4 implies that for a positive constant C2 we 
have that \9q\, \9 — 7r| > 2c2p 1//4 . This and \a + 2cos0| < 8p 3 ^ implies that we also have 
\9\, \9 — 7r| > C2P 1 / 4 provided that p is sufficiently small. On the other hand since 

2cos0 + a = -2(sinr)(0-0 Q ), 

for some r between 9 and 9$ , we deduce that for some positive constant c\ , we have that 
\9 — 9 \ < Cip 1 / 2 . Also, there exist positive constants c 3 and c 4 such that if 9 is close to 9 r , 
then 

|l — \x + y\ I > cs\a 2 + 2a cos9\ > c^^fp\9 — 9 r \, 

for r = or 1. From this we learn that if \x\ > p 1//4 , then the term is bounded above 

by 

/ [1108(^10-001)1+1]^ + —/ [|log( C4 Vp|0-0i|)i + l]^- 



2tt 



As a result, |I/(x)| < C^fp\ log p| . This in turn implies that 

#2221 < Cyp|logp| (2.36) 
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We now turn to -^2222- For this, observe that the support of the function L 2 is contained 
in the set of points x for which |x| < p 1//4 . Note that if \a\ < p 1 ^ 4 , then we can find a 
positive constant c 5 such that \9 — n/2\ < c^p 1 ^ , \9i — 3n/2\ < c 5 p 1 ^ , where 9\ = 2n — 9 . 
Forthermore, we can find a positive constant cq such that if 9 G (0, it) , then 

J a 2 + 2a cos 9\ = |2a(cos#o — cos#)| = 2 

The same is true if 9 G [n, 2tt] (use 9\ in place of 9 in ()2.H7j) .) As a result, 

|l — \x + y\\ > c 3 \a 2 + 2a cos 9\ > c 3 c e \a(9 — 9 )\. 

From this we learn that indeed 

L 2 (x) < C|log|x||l{|x| < p 1/4 }. 

To bound -R2222, let us apply Lemma 2.2 one more time to write -R2222 < -R22221 + -R22222, 
where 

-R22221 = C\ loge| 2 E N ^2 (r(xi - Xj) + r(xj - xj + zj^mirrij 

and 

-R22222 = I log e I 2 E N J (r(xi - Xj) + f (xi - Xj + z]^m i mj{d(m i ) + d(nij))dt, 

where the function r : R 2 — > [0, 00) is very similar to the function H that appeared in the 
fifth part of Lemma \2. 41 ( except that p in the definition H is replaced with p 1//4 ), and where 
the function r : R 2 — ■> [0, 00) is given by 

f(x)= / \{\x + y\< p l/4 }\\og\x + y\\dS(y). 
J\v\=i 

As in the fifth part of Lemma 2.4 we show 

iW <Cp|logp|. (2.38) 

In just the same way that we bounded J in the subsection 3.3.1, we can readily show that 
T(x) < Cyfp\ logp| . This in turn implies 

iW <Cp 1/4 | log p|. (2.39) 



9 + 9 n 



9-9 n 



a sin ■ 



■ sin ■ 



> c 6 |a( 



(2.37) 
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Putting all the pieces together we learn from (l2~3Ul) - (l2~5TJ) and (P?33j) - (l2~3Hj) that 



Ejv|^ 2 (*)|rf* < C 



l -4 + p+ \z\ + p 1/4 |lo, 

p 2 1 



for p < 1 . By making the choice p = | z\ 3 , we find that 



/ Ejv|/r 2 (*)|d*< log |^| 

Jo 



2.3.3 The case of G z (l) - G (l) 
We now estimate the term 



/ E N \G Z (1)-G (l)\(t)dt. 
Jo 



To ease the notation, we do not display the dependence of J and J on the variable t. Note 
that 



where 



nT 

/ E N \G z (l)-G (l)\(t)dt<y2D u 
Jo i=i 

n I dt^2 a(m k ,mi)V e (x k - xi)\J(x k ,m k )\ 
^° k,iei q 

e^^T^m,) 



(2.40) 



iel 



tfiXk —Xi+Z) - U (x k - Xi 



each of the other seven terms on the right-hand-side of ()2.40j) differing from Di only in an 
inessential way. Given this, the estimates involved for each of the eight cases are in essence 
identical, and we examine only the case of D\. We write Di = D 1 + D 2 , where we have 
decomposed the inner i -indexed sum according to the respective index sets 

{i G I g ,i 7^ k, I, \x k — Xi\ > p} and {i G I g ,i ^ k, I, \x k — Xi\ < p}. 

By the second part of Lemma I2.4| we have that 

C\z\ \ log el 



D 1 < 



P 



-E N / dt a(m k ,mi)V e (x k - xi), 
^° k,iei a 



where we have also used the fact that the test functions J and J are each supported on 
the set of particles of respective masses Mi and M 2 , and the fact that the total number of 



29 



particles living at any given time is bounded above by Z\ loge . From the bound on the 
collision that is provided by Lemma \2.1\ it follows that 



Note that D 2 is bounded above by 



< 



CE N dt a(m k ,mi)V e (x h - xi)l{m k = Mi} | loge | ^ 

^° k,lEl q iGlg 

t{\xi - x k \ < p}l{rrii = M 2 } u e (x k -Xi + z) - u e (x k - x { ) 1 1 J{x h m») | 
CE N dt ^ a ( m k,mi)V e (x k - xi)t{m k = M x }| loge| ^ 

- x k \ < p}t{mi = M 2 } \u e (x k -Xi + z)\ + \u e (x k - Xi)\ \ I J(x h m. 



Note that the last expectation is bounded by Lemma [2~3l because, by our assumption on a, 
we can find 7 such that a < 7, with 7 satisfying the assumption of Lemma 2.3. The upper 
bound provided by this Lemma in this particular application is computed in the last part of 
Lemma El We find that D 2 < C{p + |*|) log (p + \z\) . 
Combining these estimates yields 

D 3 < D^D 2 < C^ + C(p+ \z\)\ log {p+\z\)\. 



Making the choice p = \z\% leads to the inequality D 3 < \z\? | log \z\\ . Since each of the 
cases of [Di : i G {1, . . . , 8}} may be treated by a nearly verbatim proof, we deduce that 



T 

EJG 2 (1) - G (l)\(t)dt < C\z\hog\z\. 



2.3.4 The case of G z {2) 
Recall that 

G z (2) = -| loge| 2 ^2 <y(m k ,mi)V e (x k - xi)u e (x k - xi + z)J{x k ,m k )J{x l ,m l ). 
k,iei q 

If k, I G / g satisfy V e (x k — x{) 7^ 0, then \x k — xi\ < Roe, and so 

\x k — xi + z\ > \z\ — i?oe > kl/2, 
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provided that \z\ > 2RqE. This implies that 



l-U^Xfe — X\ + z)\ <C 



\og\x k - Xi + z\ 



< c 



log \z\ 



where in the first inequality, we used the first part of Lemma f2 .41 (restated). Applying this 
bound, and using the fact that the test functions J and J have compact support, we find 
that 

/ E N \G z (2)\dt < C log \z\ | loge| E N y~] a(m k ,mi)V e (x k - x{)dt 
Jo Jo k,iei q 

whose right-hand-side is bounded above by C\ log \z\ 1 1 loge| 1 , according to Lemma 12.11 
That is, 

/ |C^(2)|dt < C| log I^HI logep 1 . 
Jo 

2.3.5 The case of E N \X Z - X \ 

We now turn to En\X z — Xq\. Assume that \z\ > Roe. Using the second part of Lemma 
2.4, we have that 

E N \X Z - X Q (0)\ < C\z\ J J ^ h Ml (x)h M2 {y)\x - y\~ l dxdy, 

where L is a bounded set that contains the support of J and J. Using our second assumption 
on the initial data h n we obtain the bound C\z\ for Eat|X 2 — Ao|(0). 

2.4 The martingale term 

This section is devoted to proving the estimate ()2.8|) . Note that 

M Z {T) =X z (q(T),T) - X z (q(0),0) - £ (jfi+^j X z (q(t),t))dt 
is a martingale which satisfies 

3 t 1 

M 2 (T) 2 j =Y,E N [ A i (q{t),t)dt, 



E 



4=1 



where Ai(q,t) and A 2 (q,t) are respectively set equal to 
2|loge|~ 4 rf ( M i) 

idI q ,mi=Mi 



V x , u ^ Xi ~ x i + z ) J ( Xi ' Mi,t)J(xj, M 2 , *) 

jeI q ,mj=M 2 
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and 



2|loge| 4 rf ( M2 

ieI q ,mi=M2 



V,, u e (x j -x i + z)J(x j ,M 1 ,t)J(x i ,M 2 ,t) 

jeI q ,m j =M 1 



while A 3 is given by 



loge| 4 <y{mi, rrij)V e (xi — Xj) 



(2.41) 



i y-v r — 7?^ — ^ - x fc + z)J(xi,mi + m j )J(x k ,m k ) 



+ 



mj + m 



+ 



+ 



rrii + m 

m j 
mi + m 



-u e (x k -Xi + z)J(x k ,m k )J(x i) m i + mj) 
-u e (xj — x k + z)J(xj,m,i + m,j)J(x k , m k ) 
-u e (x k - Xj + z)J(x k) m k )J(xj ) m i + rrij) 



u e (xi - x k + z)J(x i) m i )J(x k) m k ) - u e (x k - x { + z)J(x k ,m k )J(x i ,m i ) 
u e (xj — x k + z)J(xj,mj)J(x k ,m k ) — u e (x k — Xj + z)J(x k ,m k )J(xj,mj) 



u 



(xi - Xj + z)J(xi, mi)J(xj,mj) 



Recall that, by our convention, we do not display the dependence of J and J on the t- 
variable. To bound these terms, we require two variants of Lemma 2.3 : 

Lemma 2.5 There exists a collection of constants C : N 2 — > (0, oo) such that, for any 
continuous functions t, v, a±, a 2 , a 3 : IR 2 — > [0, oo) and any z G M 2 , 



E 



v 



/ (it ^ m,j)t(xi - Xj + z)v(xi — x k + z) 

J ° i,j,kel q (t) 

a 1 (x i )a 2 (x j )a 3 (x k )l{m i = n x ,m k = n 3 } 

,ri2,n3 "Eji ^k)i 

"2 ij',fceig(o) 
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where A e nin2nz : R 6 -> [0, oo) zs given by 

K^x^x,) = c (6H ni ,n 2 ) £ (^^- + + ^l^ 1 ) 

t(z' - y + z)u(z' - y' + z)a 1 (z')a 2 (y)a 3 (y')dz'dydy, 

with 7 as in Lemma 2.3. 

Lemma 2.6 There exists a collection of constants C : N 3 — > [0, oo) suc/i inai, /or ant/ 
zG R" 1 , any continuous functions v, w : R 2 — > [0, oo) and another (ai, a 2 , 03) : R 6 — > [0, 00) ? 

dt ^ l(n i ,n j )V e (x i -x j )v(x i -x k + z)w(x i -xi + z)a 1 (x i )a 2 (x k )a 3 (xi) 
l{ m . = m,m fc = n 3 ,m, = n 3 } 

< Cni,n 3 ,n 3 | l°g e | ^ ^jy ^ B € m . m . mkmi {Xi,Xj,X k ,Xi) 
n 2 i,j,k,lel g (p) 

t{mi < m,mj < n 2 ,m k < n 3 ,m l < n 3 }, 
where ££ 1)n2>n3)n3 : R 8 -> [0, 00) is given by 

B n u n 2 ,n 3 ,n 3 ( x l, x 2, %3, X3) 

f / I "12 I /l2 I 12 I / 1 2 X — 3 

° Jm? V ^("i) d(n 2 ) d(n 3 ) d(n 3 ) / 

7(ni, n 2 )K(5 - z')v{z -y + z)w(z' - y' + z)a 1 (z)a 2 (y)a 3 (y')dzdz'dydy', 

with the function 7 : N 2 — > (0, 00) satisfying 

ra 2 7(ni, n 2 + n 3 ) max jl, [ ~~^~y^ ] } ^ (^2 + ^3)7(^1, ^2)- 

The proof of Lemma 2.5 is identical to that of Lemma 2.3. The proof of Lemma 2.6 is 
very similar to the proof of Lemma 2.3 and is omitted. 

We now bound the three terms. Of the first two, we treat only Ai , the other being 
bounded by an identical argument. By multiplying out the brackets appearing in the defi- 
nition of A\ , we obtain that this quantity is bounded above by 

C I log e J 4 J Vn e | (xi — Xj + z) | Vn e | (x, — x k + z) J 2 {xi, m,) 

i,j,kel q 

I J(xj : mj) || J(x k: m k ) 1 1 {rrii = M 1 ,m j = m k = M 2 ) (2.42) 
+ C|loge| 4 ^2 \u e {xi - Xj + z)\\u e (xi - x k + z)\\VJ{xi,mi)\ 2 

i,j,kel q 

I J(x j ,m j )\\ J(x k ,m k )\l{m i = M 1 ,m j = m k = M 2 ) . (2.43) 
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Let us assume that z = because this will not affect our arguments. We are required to 
bound the quantity appearing in the statement of Lemma for each of the following cases: 



(t,v, ai ,a 2 ,a 3 ) G j(|W|,|W|, J 2 , \J\,\J\) , (u e ,u e ,\VJ\ 2 , \j\, 



(2.44) 



Recall that each of the test functions J, J, and their gradients, is assumed to be uniformly 
bounded with compact support. To each of the two cases, Lemma 12.51 applies. For either 
of them, the right-hand-side of the inequality in Lemma 12.51 may be written as a finite sum 
of the expectations appearing there, with the sum being taken over triples of given masses 
ni,ri2 and n 3 . Such an expectation is bounded above by 



C loge 



-4 



l\2 



F2 - y\ 

d(n 2 ) 



+ 



E// 

n 2 m J Jk 

t{z — y)v{z — y')h ni {x\)h n2 {x 2 )h nz (x 3 )dz' dydy dx\dx 2 dx 3 , 



F3 ~ V_ 

d(n 3 ) 



/|2 



(2.45) 



where K = {x : \x\ < £} C ]R 2 is chosen to contain the support of J and J. As in Section 
3.4 of [2], we can use our bounds in the first two parts of Lemma 2.4 and repeat the proof 
of the eighth part of Lemma 2.4 to obtain 



E 



/ [A l (q(t))+A 2 (q(t))]dt<C\loge\- 1 . 
Jo 

We must treat the third term, A 3 . An application of the inequality 

(ax + ... + a n ^j < n(a\ + . . . + a 2 n 
to the bound on A 3 provided in (j2.41j) implies that 



(2.46) 



A 3 (q) <9|loge| 
where Yy is given by 



^ a(mi,mj)V e (xi - Xj) 



r 8 



E Em +n 



n=l k£l„ 



(2.47) 



rrii 



rrti + rrij 



u e (x,i - x k + z)J(xi, rrii + m j )J(x k ,m k ), 



and where {Yi : i G {2, . . . , 8}} denote the other seven expressions in (j2.41j) that appear in 
a sum over k G I q , while Yg denotes the last term in (|2.41|) that does not appear in this 
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sum. There are nine cases to consider. The first eight are practically identical, and we treat 
only the fifth. Note that 

| log e| ~ 4 ^ a(mi, m^V^Xj - xj) ( ^ Y B ) 

i,jelq k&Iq 

= C|logep 4 ^ a(mi,mj)V e (xi - Xj) 

u e (xi - x k + z)u e (xi - xi + z)J 2 (x i ,m i )J(x k ,m k )J(xi,m l ) . (2.48) 

k,l£l q 

In the sum with indices involving k, I G I q , we permit the possibility that these two may be 
equal, though they must be distinct from each of i and j (which of course must themselves 
be distinct by the overall convention). 

Note that the expression ()2.48|) appears in the statement of Lemma 12.61 provided that 
the choice 

(^v,w,ai,a2,a 3 S j = (\u e \,\u e \, J 2 ,\j\,\j\^J 

is made. Again we set z = because this does not affect the estimates. Given that the 
support of each of the functions ai, ci2, a 3 : M 6 — > [0, oo) is bounded, we must bound 

[ [ ( \ Xl ~z\ 2 | \x 2 -z'\ 2 | \x 3 -y\ 2 | \x 3 -y'\ 2 y 3 v f z - z' 
Jl*\ ^(m) d(n 2 ) d(n 3 ) d(n 3 ) J \ e 

\u e {z — y)\\u e (z — y')\h ni (xi)h n2 (x2)h n3 (x 3 )h n3 (x 3 )dzdz' dydy' dx\dx2dx 3 dx 3 , 

for a compact set L. This expression is bounded above by 



L 3 



/ z — z f \ 

Vy ) ue (z ~ y) ut {. z ' — y')dzdz dydy' 



n 2 £N 



X\ — z\ 2 \x 2 — z'\ 2 \x 3 — y\ 2 \x 3 — y'\ 2x 



d{ni) d(n 2 ) d(n 3 ) d{n 3 ) 

h m {xi)K 2 (x 2 )h n3 (x 3 )h n3 (x 3 )dx 1 dx 2 dx 3 dx 3 , 

which is less than 

C J v(^-^ju e (z-y)u e (z' - y')dzdz' dydy' . 

The proof of this follows the proof of the eighth part of Lemma 2.4; we use the elementary 
inequality abed < (a 2 + b 2 + c 2 + d 2 ) 2 and the fact that the function 

k{x) = Y,d{nf" j h n (y)\x - y\^ 2 dy 
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L 



is locally bounded. Noting that the bound < | log|x|| implies that 

u £ (z-y)dy 
is bounded above by a constant, we find that 

J V^—^)u e {z - y)u e {z' - y')dzdz'dydy' <C J V^—-^jdzdz'. 
This is at most Ce 2 . Applying Lemma \2. 61 we find that the contribution to 

E N [ A 3 (q(t))dt 



o 



arising from the fifth term in (j2.47J) is at most 

Ce-^V^Ilogel" 1 . 
We now treat the ninth term, as they are classified in (J2.47j) . It takes the form 

| loge| 4 ^2 a { m u m j) v e(xi 



u e [xi — Xj + z) 2 J(xi, mj) 2 J [xj^nij) 2 . 



This is bounded above by 

C|loge| 3 a(mi,mj)V e [Xi - Xj), 



because u £ < C\ loge in the support of V £ by the first part of Lemma 2.4. The expected 
value of the integral on the interval of time [0,T] of this last expression is bounded above 
by 

3 f T I 

C|loge| Ejy / dt a(mj, mj)V £ (xj — Xj) < C \ loge| , 
where we used Lemma 12.11 for the last inequality. This completes the proof of (|2.8|) . 
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2.5 Using the estimates 

The inequalities (|2.fij) . (|2.7jl and (|2.8jl imply that, for large T, 



lim limsupEAT 



H n (t)dt + / H 13 (t)dt 

i! >/0 



di = 0. 



That is, 



lim lim sup I log el Ejv 



a(mi,mj)J(xi,mi)J(xj,mj) 



V e (xi — Xj + z) — V e (xi — Xj) + V e (xi — Xj + z)u e (xi — Xj + z) 



0. 



(Recall that we simply write J{xi 1 m i ) and J(xi,rrii) for J{xi,rrii,t) and J(xj, m i? t) .) This 
implies that 



loge 


7 




</ 


loge 


7 




Jo 



o „• 



a(mi, rrij)V e (xi — Xj)J(xi, ra i )J{xj 1 rrij)dt 



*j'G/ 9 (t) 



where Err\ satisfies 



a(rrii, rrij)V e (xi — a;^ + z)J(xi, mi)J(xj, rrij) 
|^1 + | log e| 1, u e (xj — Xj + z) dt + Erri(e,z), 



lim lim sup Ejv I £?rri (e, z) I = 0. 



(2.49) 
(2.50) 



By Theorem 3.2, the expression 1 + | loge| 1 u e mi m2 (a) is uniformly close to 

r(mi, m 2 ) 



1 



27r + r(mi, m 2 ) 



for a satisfying K e (a) ^ and r(m 1 ,m 2 ) = a(mi,m2)/(d(mi) + d(m 2 )) . Recalling from 
([□1 that 

Q=|loge| (^(m^mj^K^i - Xj)J(x i ,m i )J(x j ,m j ) 

and writing 



Q(z) = | log e | 1 (3(mi,mj)V £ (xi Xj + z)J(x i} m i )J(xj,mj), 



(2.51) 
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it follows from (|2.49|) and Theorem 3.2 that 



where Err 2 satisfies 



Q(t)dt = / Q(z)(t)dt + Err 2 (e,z), 
Jo 



lim lim sup Ejv I Err 2 (e, z) I = 0. 
l*l-o 40 1 n 



(2.52) 



From this, it is not hard to deduce that 
Q{z 2 - z x ) = 



loge| 2 P( m i, mj)V £ (xi - xj + z 2 - z\) 

J(xi - z 1 ,m i )J(x j - z 2 , rrij) + Err(e, z 1: z 2 ), 



(2.53) 



where 



E N \Err(e, Zl , z 2 )\ < C(\ Zl \ + \z 2 \). 
(See Section 3.5 of [2 .) By (Z5Z> and (g33D , 



T 



Q(t)dt 

loge|~ 2 / dt P(mi,mj)J(xi — Zi,rrii)J(xj — z 2 ,rrij) 

X X ^ (( ^ " Zl) ~ (Xj ~ Z ^ 5 ~ 2r l(j) 5 ~ 2 v(j)dzidz 2 + Err 3 (e,5) 

log£|^ 2 / dt dujiduj 2 \ V £ {uj\ — LU 2 )(3(mi 1 rrij)J(iUi,mi)J(uj 2 ,mj) 
Jo Jr m 



r 2 r/ 



Xi — UJ\ \ 2 / -Ej — ^1 



5 



6 



T 

dt 

o Jm 2 




+ Err 3 (e, 5) 
duo x du 2 v E {u x - U2)/3(Mi, M 2 )J{u x , Mi)J(w 2 , M 2 ) 



loge| 1 ^ 5 2 ??( f . 

iaI q :m,i=Mi 

logep 1 ^ 



j£lq-,mj=M2 



+ Err 3 (e, 5). 



where i?rr 3 satisfies 



lim lim sup E^v I Err 3 (e, 6) I = 0, 

810 £|0 
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and where in the last equality, we made use of the fact that the test functions J and J take 
non-zero values only on particles of a given mass, respectively Mi and M 2 . Thus, 

Q(t)dt 

cT 



dt / du(3(M 1 ,M 2 )J(u;,M 1 )J(uj,M 2 ) 

JR 2 



loge |- E *-,(*^) 



i&I q :m,i=Mi 

_l_ Err(e, 5) 



where Err = Err 5 + err with the function err = 0(e5 5 ) also satisfies 

lim lim sup K N I err (e, S) I =0. 

<5|0 ei o 



This completes the proof of Proposition I. 



3 Potential theory 

The purpose of this section is twofold. Firstly, we show the existence of the function u e that 
satisfies (3.1). Secondly we evaluate the limit of u e \ loge| _1 in the support of V s , as e — > 0. 
This limit was used in the evaluation of (3 in Section 2.5. We start with the statements 
of the main results of this section. Let V : M. d — > R be a continuous function of compact 

support with V > and / ^(x)^ = 1. We also write K for the topological closure of 

C/o where 

U = {x: V{x) ^ 0}. (3.1) 

Given a measure /x, let us define 



Gn(x) = J log \x - y\n{dy). 



When the measure \x is absolutely continuous with respect to the Lebesgue measure with a 
density g, we simply write Qg for Q\i. 

Theorem 3.1 There exists a number 70 > such that for every 7 e (0,70) and a G M, 
t/iere exzsfe a unique function u E C 2 (M 2 ) snc/i £/ia£ = 0(| log|x||) as \x\ — > 00 and 

M = 7 £(( M + a )\/). (3.2) 

Moreover Z = 7 f(u + a)Vdx 7^ and (n + a)Z _1 > . 
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Recall that we are searching for a function u £ such that 

Au 6 = r(n, m) [V e u e + V £ ] , 
where r = r(n, m) = a(n, m)/(d(n) + d{m)) . For this it suffices to have 

u £ = Gfi £ , (3.3) 
where fi £ (dx) = ^r(V e u £ + V £ )dx. This can be rewritten as 

u £ = ^r £ Q(y £ [«- + | loge|]) , 

where r £ = r| logej" 1 . Evidently we can apply Theorem 3.1 to deduce the existence of the 
function u £ for sufficiently small e. 

Our next theorem was used in the previous section for the evaluation of (3. 



Theorem 3.2 For every positive k , 



lim sup 

\x\<k 



u £ (ex) \ log el 1 + 



T 



2n + t 



Proof of Theorem 13.11 

Step 1. Let J be a bounded continuous function with J > and 



J(x) (log \x\) 2 dx < oo. 



Define 



H = < u : u is measurable and / u 2 {x)J{x)dx < oo > . (3.4) 
I Jm 2 ' ) 

We then define JF : 7i — >• TC by ^(u) = Q{uV). Observe that 7i is a Hilbert space with 

respect to the inner product 



(u, v ) = / u(x)v(x)J(x)dx. 
Jr 2 

Let us verify that T is a bounded operator. To see this, write 

T(x) = / I log \x - y\\V(y)dy. (3.5) 
Jr 2 
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When \x\ is sufficiently large, we have that T(x) < log(2|rr |) because < log | x — y\ < 
log(2|x|) whenever V(y) 7^ 0. Otherwise we have 



r(z) < c / I log [a; -y\\dy < c 2 , 

J\x-y\<c\ 

for constants c ,Ci and c 2 . As a result 

r(x) < c + iog + \x\ 

for a constant c. Also, we may use Holder's inequality to assert 



(F(u)(x)f < 



Y{x) / I log |a; - y\\V{y)\u{y)\ 
Jr 2 



dy 
T(x) 



< T(x) I log|s - y\\V{y)u 2 (y)dy. 



From this we deduce 

{T{u){x)fJ{x)dx< I V(y)u\y) 



T(x) I log \x — y\ I J(x)dx 



dy. 



If V{y) 7^ then |y| < _R for a suitable _R . Define 

I(y) = / r(x) I log \x — y\ I J(x)dx, 



Note 



+ / r(a;) I log \x — y\ | J(x)dx 

|a;|<2_R J\x\>2R 



(3.6) 



(3.7) 



< Ci / I log |s — j/| \dx + ci / I log |s| | 2 J(x)dx, 

J\x-y\<3R J\x\>2R 

where for the second line we have used (J3.6|) . From this and our assumption on J we deduce 

that sup I(y) < 00. As a result, 

\y\<Ro 



{F{u){x)yj(x)dx< Cl V{y)u\y)dy<c 2 u\y)J{y)dy 
r 2 Jr 2 Jr 2 

because V is of compact support and J > 0. This shows the boundedness of the operator 
T -.H—tH. 

Step 2. Since the operator JF is bounded, the equation 

(id - iF){u) = g 
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has a solution, where g(x) = — , ~faT(x) with T as in ()3.5|) and id denotes the identity 
transformation. Note that our assumption on T implies that V <EH because of (13. 6|) . So far 
we have shown the existence of a unique solution u G Ti of u — rjF(w) = g . From this and 
the Holder continuity of V we can readily show that in fact u G C 2 and that u is a classical 
solution of 

Am = 2-k 1 (u + o)V. (3.8) 

(See for example Section 4.2 of 0.) 

Step 3. In this step we verify Z^O. Observe that u = Q^i for a measure « with a bounded 
support. From this we can readily deduce 



u (x 



) = /i(M 2 )log|a;|+0(|a;|- 1 ), (3.9) 



Vu(x) = fi (M 2 ) + O (\x\- 2 ) . (3.10) 
We now choose R> Rq and use Am = 2ttj(u + a)V to write 

(u + a)Audx = 27T7 / V(u + a) 2 dx. 
\x\<R J\x\<R 

After an integration by parts we obtain 

— / \Vu\ 2 dx + I (u + a)Vu ■ ndS = 2^ / V{u + a) 2 dx, 

J\x\<R J\x\=R, J\x\<R 

X 

where n = - — - is the normal vector and dS is the Lebesgue measure on \x\ = R. Now if 

\x\ 

Z = fi(R 2 ) = 0, then we can use (l3~H- (l3~TTH) to deduce that 

I {u + a)S7u-ndS = 0(R- 1 ). 

J\x\=R 



As a result, 

- / \Vu\ 2 dx = 2vr7 / V{u + a) 2 dx. 

From this we deduce that / \Wu\ 2 dx = i (u + a) 2 Vdx = . This in turn implies that 

Jr 2 Jr 2 

u = 0. But this contradicts u = Q(V(u + a)). Hence we can not have Z = 0. 
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Step 4. It remains to show that (u + a)Z~ 1 > 0. We only establish this when Z > 0. The 
case Z < can be treated likewise. First take a smooth function (p$ : R — > (— oo, 0] such 
that <// 5 > and 

r > —a, 
a + r r < —a — 5. 



We then have 



/ (fs(u)Vu.ndS — / <p' s (u)\Vu\ 2 dx = 2^ / ^(^(u + ajFrfj; (3-11) 

J\x\=R J\x\<R J\x\<R 



by an integration by parts. (Here n = x/\x\.) If Z > 0, then we can use ()3.9j) to assert that 
u(x) > and <ps(u(x)) = whenever |x| = R and R is sufficiently large. Since the left-hand 
side of (j3.11J) is negative for such large R and (u + a)(p${u) > we deduce 

ip' 8 (u)\Vu\ 2 dx = I V(u + a)(ps(u)dx = 0. 
Jr 2 

We now send 5 — ► to deduce 

0= / \Vu\ 2 \{u + a < 0)dx = / V(u + a) 2 l(u + a < 0)dx. 
Jr 2 Jr 2 

As a result, on the set A = {x : a + -u(x) < 0} we have Vu = 0. Hence u is constant on 
each component B ol A. But this constant can only be —a because on the boundary of 
A we have u + a = 0. This is impossible unless A is empty and we deduce that u > —a 
everywhere. □ 
We now turn to the proof of Theorem 3.2. We first state and prove a lemma. Let us 
write A e for //(R 2 ) where // was defined right after (3.3). 

Lemma 3.1 We have that A e > for small e. Moreover 

limsupA e < 1. 

Proof Let us write u e for {u e + jlogeDAj 1 and fi s for A^ 1 //. By Theorem 3.1 we 
have that fi e is a probability measure and u e > . Note that the support of the probability 
measure fi e is the set eK . Moreover u e is harmonic off eK and 

if = + i^fl. (3.12) 
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By a well-known theorem in potential theory we have 

|loge| 



eCap(K ) = Cap(eK ) > exp 



A, 



where Cap denotes the logarithmic capacity (See for example Theorem 9.8 of As a 
result 

logCap(A-o) x> 1 
|loge| ~~ A e ' 

From this, we can readily deduce the claims of the Lemma. □ 
Proof of Theorem 3.2 It suffices to show that for every positive k, 



lim sup \u £ (ex)\ \ogs\~ l + Aj = 0, (3.13) 

and 

limA e = — - — . (3.14) 

e-o 2tt + t y 1 

Recall that by Theorem 3.1 and Lemma 3.1 the expression u £ (ey)\ loge| _1 + 1 is nonneg- 

ative Also recall that Rq is defined so that the ball Br o (0) contains the support of V. Let 
us write 4 for the maximum of u £ (ex) over the ball 5^(0). We then have 



u%ex) = 1- 



J log \ex — ey\ (u 6 (ey)\ log e| 1 + l) V{y)dy (3.15) 
A e loge + ^- J \og\x-y\ (u e {ey)\\oge\~ l + l) V{y)dy 



< A £ loge + ^- (4|loge| : + l) / \og\x -y\V(y)dy. 

ZlX J\x-v\>\ 



\x-y\>\ 

Hence for every x with \x\ < k we have 

u £ (ex) < A £ logs + c (4| loge]" 1 + l) , 
where c is a constant that depends on k only. By choosing k = Rq we deduce that 

(l - c\ loge| _1 ) 4 < A e loge + c. 



Hence 



4 = max u £ (ex) < 2A £ \oge + 2c (3.16) 

|x|<Ro 



for sufficiently small e. This in turn implies 

4|loge|" 1 + 1 < 2, (3.17) 
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for small e. Moreover, by the second equality in (|3.15|) . 

u £ {ex)\ \oge\- 1 + A £ = — f , / log \x - y\ [u E {ey)\ \oge\~ 1 + l) V{y)dy = X X + X 2 , 

2ir\ loge| J v ' 

where 

X 1 = — - j \og\x-y\ (^{ey^loge^ 1 + l)V(y)dy, 

27r| loge| J| x _,,|<i 

X 2 = — — - I \og\x-y\ (u £ (ey)\ \oge\~ 1 + l) V{y)dy. 

27r| Joge| y[ x -v[>a 

Since the expression u £ (ey) \ logej^ 1 + 1 is nonnegative and bounded above for y in the ball 
-Bf?o(0), we deduce that both X\ and X 2 converge to in low e limit. This completes the 
proof of (3.15). 

We now turn to the proof (3.16). By the definition of A e and (3.15), 
A £ = (u%ey)\loge\- 1 + l)V(y)dy 

= - A e ) J V(y)dy + o(l) = ^-(1 - A.) + o(l). 

This immediately implies (3.16). □ 
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